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Abstract 

We explore the possible interpretation of the large-Q 2 anomaly reported by the HI 
and ZEUS collaborations in terms of stop squark production off a strange quark in the 
proton via an R- violating interaction. This "strange stop" interpretation is constrained 
by LEP measurements of the Z — > e + e~ decay rate in addition to constraints from the 
electroweak p parameter and CDF and DO searches for first-generation leptoquarks. We 
investigate the interplay between these constraints, taking full account of stop mixing 
effects. We find that if < 200 GeV only relatively small domains of the chargino and 
neutralino parameters are consistent with these constraints, and explore the extent to 
which this scenario may be probed further by searches for contact interactions at LEP 2 
and experiments with e~ and polarized beams at HERA. 
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The surprisingly large number of high-x and -Q 2 deep-inelastic e + p scattering events at 
HERA, reported recently by the HI and ZEUS collaborations 0, is currently exciting con- 
siderable theoretical speculation |3j, H, II, El- ft ma y wen be that the excess reported is just a 
statistical fluctuation, and this may become apparent from the additional data being obtained 
during the current year's run. However, unless and until this happens, it is the responsibility 
of theorists to pursue different possible interpretations of the excess, assuming it to be real, 
evaluate and combine the constraints imposed by other data, and propose additional signatures 
that may serve to discriminate between the surviving rival interpretations. Effective contact 
interactions from some energy scale beyond that of the Standard Model do not seem to match 
well the combined Q 2 distribution reported by the HI and ZEUS collaborations. Interest has 
been piqued by the tendency of the HI data to peak in an x range corresponding to the produc- 
tion of a narrow resonance with a mass around 200 GeV Jl| and leptoquark quantum numbers, 
which is apparently not contradicted by the ZEUS data 0, [7]]. 

If one chooses to speculate along this line, upper limits on the production cross section 
at the Fermilab Tevatron collider || [9| strongly suggest that any such leptoquark state must 
have spin 0, and hint that its branching ratio B(e + q) into the observed final state should be 
significantly less than unity [|J. The latter is difficult, though not impossible, to arrange in a 
model without supersymmetry, but is natural if the "leptoquark" is actually a squark produced 



via an interaction that violates R parity |L(| [Ll], |12|]: \' i - k L i QjD k l , where L h Qj denote left 



handed lepton and quark doublets, and D c k denotes left-handed charge-1/3 antiquark singlets. 
In such a supersymmetric scenario, there may be significant competition between the branching 
ratio B for _R-violating decays via the A' ljfc production coupling and i?-conserving decay modes 
such as XQ followed by x ~^ ^QQ decay [||]. Within this general supersymmetric framework, 
three possible interpretations seem to emerge: e^dR — > cl, e^dn — > i and e^SR — > i. 

In this paper, we focus on the latter "strange stop" interpretation, which requires a relatively 
large R- violating coupling A' 132 ~ 0.3/(cos(/>a/S), where is an angle describing mixing between 
the tz,R Q, and possibly fine tuning of the chargino spectrum so as to have a significant, but 
not too dominant, branching ratio into x + b 0- As was pointed out in 0, this scenario must 
contend with LEP 2 constraints on virtual t exchanges from the OPAL collaboration [13], as 
well as with upper limits from precision electroweak physics on the possible contribution of 
the t,b sector to the p parameter. In addition, as was pointed out in |TJ[], the measurement 
of the Z — > e + e~ decay rate provides an upper bound on the possible magnitude of any X' 13k 
coupling, which is of particular interest because it is independent of the stop mixing angle 
and decay branching ratio. As we shall see, new data strengthen significantly the constraint 
found in the previous analysis, providing stimulus for the re-evaluation in this paper of all the 
phenomenological constraints on the "strange stop" interpretation of the HERA data, taking 
into account tL,R mixing effects 0. We do not comment on possible stringent limits on im- 
parity violating interactions from cosmological considerations [|15[], as they can be avoided in 
various schemes, such as baryogenesis at the electroweak scale |16|] . We analyze the restrictions 
imposed by the Z — > e + e~ and other constraints in the /i, M2 plane that characterizes chargino 
and neutralino properties 0, and hence stop decay modes. If the stop mass is 200 GeV or below, 

2 This must be significant, because of the electroweak p parameter 

3 We also comment, where appropriate, on implications for the "down stop" interpretation. 

4 We assume universality at the supersymmetric GUT scale for the input soft supersymmetry-breaking gaug- 
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there is a restricted region of this plane where A' 132 is sufficiently small to satisfy the Z — > e + e~ 
constraint and B(e + s) is sufficiently small to avoid the Tevatron production constraint. The 
latter imposes no constraint if the stop mass exceeds 210 GeV, in which case a significantly larger 
region of the fi, M2 plane is allowed. We find that searches for effective contact interactions 
at LEP 2 may have a limited chance to probe further the "strange stop" scenario, and also 
comment on the possibilities for e~ or polarized beams at HERA to cast light on this scenario. 
Any coupling of the form X' i3k would make an extra contribution to the Z partial width 
= Y[Z — > ifl^) via loop diagrams involving a t and a c4, which is proportional to for large 
m t . Since this may become quite large, precise determinations of the T/. at LEP can provide 
sensitive constraints on the possible magnitude of R violation in the top sector jl4) . There are 
also contributions to STi from diagrams with stops in the loops, but these are suppressed by 
a factor mf / vr? t as compared to the dominant contribution, where i — d,s for the couplings 
relevant to the interpretation of the HERA anomaly. These diagrams have negligible numerical 
impact on the bounds that we derive, so the effects of stop mixing on our bounds are also small. 
We note that the masses of the dk could in principle be different from that of the t assumed 
to be produced at HERA. However, the numerical values of the loop diagrams are not very 
sensitive to m d - fc , which we assume for definiteness to equal m^. 

We evaluate the top-c4 loops using the full expressions for their contribution 5T [t given in 
Ref. ||14|| . It is clear that the dominant contributions come from diagrams with a top quark 



and a d or s squark in the A' 131 and A' 132 cases, respectively. The non-observation of flavour- 
changing decays such as p — > imposes upper limits on the possibility that the lepton-number 
violation induced by A' couplings can be present simultaneously in both e(i = 1) and p{i = 2) 
channels: A' 132 A 232 < 0.015 for mi ~ 200 GeV |17| . This implies that these couplings cannot 
affect significantly both r e and Q, so we improve the bounds on the A' 13fc couplings presented 
in Ref. |14]|, by considering the ratio T e /T ^ []. We have used the following experimental values 



for the leptonic widths 0: T e = 83.94 ±0.15 MeV and = 83.79 ±0.21 MeV. The resulting 
bounds on A' 13fe are shown in Fig. 1. We see that the bound on A' varies from about 0.52 to 
about 0.6 as the squark mass changes from 180 to 220 GeV. We emphasize that this bound is 
independent of stop mixing and the t — > e + q branching ratio B, making it a valuable constraint 
on "strange stop" interpretation of the HERA anomaly. 

The mass splitting in the stop-sbottom sector results in a contribution to the electroweak p 
parameter, which is constrained by precision Z measurements from LEP and the SLD, and by 
mw/mz measurements from the Tevatron and elsewhere. These and the Tevatron measurement 
of mt provide strong upper limits on any contribution to p from physics beyond the Standard 
Model, in particular from the stop-sbottom sector. In the absence of stop mixing the bi 
mass is given by m~ bL = \J^i 2 L ~ cos2/3M^ / — mf ± mf, and is about 100 to 130 GeV, which 
would result in a contribution Ap larger than that allowed by the data at the 95 % C.L. 0. 
We now evaluate the amount of mixing that is required to avoid any contradiction with the 
experimental upper limit on Ap. 

We neglect the bottom quark mass, which also implies that we may neglect mixing in the 



ino masses Mi j2 . 

5 This also means that B(p + s) << B(e + s) in the "strange stop" scenario. 

6 This is equivalent to considering the ratio i?, e /i? M , because the corrections the hadronic width Th from loops 
with sleptons cancel when one takes the ratio. 
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sbottom sector. The mixing in the stop sector may be expressed in terms of an angle (j) given 
by 

m~ — m\ — m? — m 2 , / cos2/? 
tan0 = -± ^ -. (1) 

m t m LR 

where mg is the mass of the physical left-handed sbottom, and we parametrize by mtrriLR the 
off-diagonal element in the stop mass-squared matrix where nii^ is a combination of the soft 
supersymmetry-breaking A parameter and the Higgs mixing parameter /i. We denote by 
the mass eigenstates that result from the mixing: 

,2 _ 1/2 , _2 , n 2 , 2 „„„na/l , 2 + „„2/ 



m *i,2 = 2^ m ~ bL + m ~ tR + 2m * + m w cos2 K' i - + 3 tan 

T{l(m- bL - m 2 iR + m 2 w cos2/3(l - |tan 2 ^)] 2 + im 2 t m 2 LR }^, (2) 



where m 2 R is the soft supersymmetry-breaking contribution to the right-handed stop mass. The 
contribution to Ap coming from the t — b mass-splitting is given in this notation by [19|, ^0 

Ap(t -b) = — sin 2 6W \ cos 2 0/ s (m£ , mf ) 

sin 2 0/ s (m| , m~ bi ) - cos 2 </>sin 2 0/ s (m| , (3) 



where 



mkn , m» 1 
— - — — In — - + - 

2 2 2 ' r> 

mf — wij mf 2 



/ s (m 2 , m 2 2 ) = , x % ln^ + ~(mf + m 2 ) (4) 



We note that the (t, 6)-loop contribution to p is positive, as is the range of p indicated by the 
precision electroweak data. 

We plot in Fig. 2 the t, b contribution to Ap as a function of the mixing angle for different 
choices of tan/? and ttilr for = 200 GeV. The parameter ttilr cannot be much larger 



than mi L , because of the problem of false vacua |2~I |. As a rule of thumb, one may impose 
^ rn LR ^ 3m^, in which case sincp < 0. For the purposes of this analysis we choose 
wlr = kmi L , with k = 1,2. The curves marked (a) and (b) in Fig. 2 are for the case k = 1 
and tan/3 = 1,5 respectively, whereas those marked (c) and (d) are for k = 2 and tan/3 = 1,5. 
Curve (e) represents the 95 % CL. upper limit after subtracting the Standard Model reference 
value for m# = 100 GeV, as is appropriate for the MSSM. We see that for k = 1 and tan/? = 1, 
values of <fi > —0.6 are excluded for m H = 100 GeV. 

We now calculate Bit — > e + s) in the context of the "strange stop" mechanism e + s R — > t 
studied in this paper, including stop mixing effects Q In the presence of mixing, the i?-violating 
t\ — > e + s decay rate becomes 

r*! = Y^ A i32 m 9 cos2 ( 5 ) 

7 In 0, the corresponding calculation of B(tL — > e + d) was presented, neglecting stop mixing and under the 
alternative assumption that the t l was produced off valence d quarks in the proton, with a much smaller value 
of the corresponding A' 131 coupling. 
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Figure 2: Contributions to the p parameter from t—b loops, for the cases (a) k = 1 and 
tan/3 = 1, (b) k = 1 and tan(5 = 5, (c) k = 2 and tan(3 = 1, (d) k = 2 and tan(3 = 5, 
compared with (e) the 95 % C.L. after subtracting the Standard Model prediction for 
m H = 100 GeV. 
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and the decay rate of t 2 is described by a similar equation, with cos0 substituted by sin0 0. To 
illustrate our analysis, we choose values of the mixing angle <fr that are consistent with the Ap 
bound shown in Fig. 2. 

The most important i?-conserving decay is presumably that into a chargino: t — > x + b. 
To understand this, we first note that the decay t — > where \ is the lightest neutralino, 
is presumably excluded by LEP 2 lower limits on m x , though a thorough analysis in the 
appropriate i?-violating framework has not yet been published. Secondly, the kinematically- 
allowed t — > x c decay is suppressed by typical loop factors, and three-body decays are also 
expected to be small J3|. Since the value of the A' 132 coupling required in the "strange stop" 
interpretation studied here (?Z 0.3) is much larger than that of the corresponding A' 131 coupling 
in the "down stop" interpretation 0.04) pursued in 0, it is easier for the R- violating decay 
mode to become competitive with the .R-conserving chargino decay mode, even when the latter 
is not kinematically suppressed. 



In the presence of mixing, the t — > x + 6 decay rate becomes [111 



h Asm 2 6 w 1 m 2 m 2 



(\Gl\ 




m %± \ 4ra&m.c± . 



+ l^l 2 ) ( 1 - 5r - -H - — hr~^ (GrG* l ) \ , (6) 



m~ I m~ 
ti / ti 



_ m b U£ 2 cos4> m t V k2 sin(f) 

Gl = — a! Gr = v kicos<f> + -t= — (7) 

y 2m w cos p y 2m w sm.p 

where X(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2yz — 2zx is the usual phase space factor, and the Vjy 
and Uki are the chargino mixing angles, determined in the usual way by fi, the ratio tan/? of 
supersymmetric Higgs vacuum expectation values, and M 2 . As compared to zero mixing, we see 
that the vertex involving changes by a factor cos0 like the i?-violating coupling. However, 
now we also have the vertex that couples the chargino to £r, which is dependent on V12, m ti 
and sin0. Note also that the branching ratios of the stop decays to charginos depend on the 
sign of 0, which is sensitive in particular to the soft supersymmetry-breaking parameter A and 
the Higgs mixing parameter /1. 

The presence of the additional vertex suggests at first sight that large stop mixing might tend 
to favour the i?-conserving decay mode relative to the i?-violating one. Indeed, this feature is 
manifest in the region of parameter space where /1 is small. However, there is also the possibility 
of a cancellation that may amplify the i?-violating branching ratio in some specific regions of 
parameter space. To see this, note that, when stop mixing is taken into account, besides the 
squared couplings of the ti,,R vertices, there also exists a term in the decay rate arising from the 
cross product of the parts of the tn and ti, vertices that have the same helicity f\ This "cross 
term" is proportional to 

m t V k2 V kl 



2 cos<t> sine 



\^2m w smP 



^This decay would be relevant if t 2 were light enough to have been produced at HERA, as proposed in |6). 
3 The additional term involving Ke(GnG} j ) gives a much smaller contribution, except in the small region 
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This is significant, and can be destructive if either (i) sine/* is negative and V\\ and V\ 2 have the 
same sign, or (ii) sin0 is positive and V\\ and V\ 2 have opposite signs. 

It is clear that the stop branching ratio B(e + q) must be non-negligible, if HI and ZEUS 
are to have seen this decay mode without being drowned by other stop decay signatures. On 
the other hand, the DO || and CDF || collaborations at the Fermilab pp collider have both 
established upper limits on aLQB(e + q) 2 , where LQ is a generic scalar particle with decays into 
e + q, which may for our purposes be taken as the t\. In particular, CDF || establishes that 
m LQ > 210 GeV if B{e + q) = 1, corresponding to B{e+q) < 1/^1 + 91n (210 GeV/m LQ ) for 
values of the leptoquark mass below 210 GeV. The CDF limit |§ corresponds, in particular, to 
B{e + q) < 0.73(0.87) for m LQ = 190(200) GeV, but provides no constraint for m LQ > 210 GeV. 

We combine the above constraints in the /x, M 2 plane as shown in Fig. 3, using the following 
procedure. We first choose values of tan/? and of </>, the latter compatible with the Ap upper 
bound shown in Fig. 2. We then use @ and (0) to calculate B(e + s), employing self-consistently 
the value of A' 132 suggested by the HERA observations: A' 132 = 0.3/(cos</>\/F). We then apply 
the Z — > e + e~ constraint of Fig. 1 - which excludes models with low B, in particular at small M 2 
and/or and the CDF constraint [f| - which excludes models with large B if < 210 GeV. 
As we see in Figs. 3(a,b), only restricted regions of the /i, M 2 plane are compatible with these 
constraints if =190 GeV, and these are not greatly expanded if one chooses mi x = 200 
GeV as in Figs. 3(c,d). On the other hand, it is clear that the entire upper part of the /i, M 2 
plane is allowed if > 210 GeV, since the CDF constraint || then imposes no restriction. 
Comparing Figs. 3(a,c) with Figs. 3(b,d), we see that the allowed regions are of similar size, 
but more symmetric between positive and negative /i, if one chooses tan/3 = 5. 

As was pointed out in [Q, one of the potentially interesting constraints on leptoquark or 
i?-violating squark interpretations of the large-Q 2 HERA anomaly could come from LEP 2 
limits on effective contributions to e + e~ — > qq due to crossed-channel exchanges. Formulae 
have been given in ||] and elsewhere for the case of single leptoquark or i?-violating squark 
exchange. However, the LEP 2 constraint is potentially most relevant for the "strange stop" 
interpretation with its relatively large coupling: A' 132 ~ O.3/(cos0v / /5), in which case stop 
mixing should be taken into account in estimating the exchange contributions. Accordingly, 
we now present for completeness the simple generalization of the standard formulae available 
previously to the case in which two mixed stops t\ y2 are exchanged: 



3A' 4 (/ lC osV + / 2 sin 4 0) 3A' 4 (/ 3 cos 2 (/>sin 2 </>) 

SAI 647TS 647TS 



3A /2 a em (/{cos 2 </> + I^sm 2 ^ 
As 



e e e d + a L a R 



d s(s - Mf) 



(9) 



where 



= TT^-- 2 ^ ln (— j (10) 

/, = i + x l m (1±Sl\ — l—hfl+a) (ii) 

X h + x i 2 \ x ii J x h + x i 2 \ X t 2 J 



n, 2 = ~ + x ^- x L^[—r Jhl ) (12) 
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Figure 3: Regions of the /x, M 2 plane that are consistent with the upper bound on A' 132 shown in 
Fig. 1 and the CDF upper limit on B(e + q) || for the i?-violating decay of ti, for (a) sin0 = —0.7, 
tan/3 = 1.0 and = 190 GeV, (b) sirup = —0.7, tan/3=5.0 and the same value of m^, (c) 
sin0 = —0.7, tan/3 = 1.0 and = 200 GeV, and (d) sin0 = —0.7, tan/3=5.0 and the same 
value of m f - . 
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with a LR = (T3 — sin 2 6Vk) / (sin 6* w cos 0^) an d 2 = ~ t" 2 "- The contributions proportional 
to Ji j2 in eq.(^) arise from the squared amplitudes of the diagrams with the R- violating vertices, 
for i\ and t 2 exchange respectively, and J3 is the interference term between the t\ and t 2 exchange 
diagrams. Finally, the terms proportional to I[ 2 are interference terms with the Standard Model 
s-channel 7 and Z exchanges, for both stop exchanges 0. 

The interference of the new diagrams with the Standard Model processes is in our case 
destructive. In the case we discussed previously [|J with A' 131 = 0.04, interference gave the 
dominant deviation from the Standard Model prediction. However, in the case studied in this 
paper with A' 132 ^ 0.3, the squares of the new diagrams are also large, and cancellations occur 
between them and the negative interference terms. Moreover, when mixing is introduced, the 
squared terms and the destructive interference ones are multiplied by different factors, in such 
a way that the cancellation is enhanced. Thus the bounds derivable from LEP 2 bounds may 
be relaxed significantly compared to the unmixed case. Indeed, for any given value of A' 132 , 
there is even a value of <ft where no bound is obtained, since the cancellation becomes exact. 
More generally, when two light stops are exchanged, the cancellation occurs over a wider range 
of mixing angle 0, and the LEP 2 bounds are relaxed. The numbers for some representative 
parameter choices are shown in the Table. For various different values of m~ t2 and <ft, the 
magnitudes of the destructive interference terms and the new squared amplitude contributions 
are shown, and we see how they conspire to leave small net contributions to the cross section 
for e + e~ — > ss. 

We have seen in the above analysis that the "strange stop" scenario may survive the available 
constraints, at the price of requiring a relatively restricted region in the /i, M 2 plane if m~ tx < 200 
GeV. It is disappointing that partial cancellations between negative interference and positive 
squared-amplitude terms may hinder attempts to pin down this scenario at LEP 2. In principle, 
the decay modes and branching ratios observable at HERA could discriminate between the 
"down" and "strange stop" scenarios. There are also possibilities for discriminating between 
these scenarios by using different beams at HERA. 

One clear distinction could be drawn by comparing stop production in e + p and e~p collisions. 
In the "down stop" case, one predicts a(e + p) » cr(e~p), since production occurs off valence 
quarks in e + p collisions, and off sea quarks in e~p collisions. One does not expect such a large 
ratio in the "strange stop" interpretation, since production occurs off sea quarks with either 
beam. However, it is not necessarily the case that the production cross sections must be exactly 
equal, since the s and s distributions may have different shapes, though their integrals must 
be equal. If the strange sea in the proton is extrinsic, e.g., generated by perturbative QCD, 
one would expect strict equality in the cross sections, but this need not be the case if there 
is an intrinsic |22|] strange component 0. Thus, we predict only that a(e + p) ~ a(e~p) in the 
"strange stop" interpretation. 

Another tool for discriminating between models could be provided by polarized beams at 



10 Note that, as well as the different coupling factors, because of fermionic Wick ordering there is a difference 
in the sign of the interference term for charge-1/3 qq final states such as the ss case considered here, as compared 
to charge-2/3 final states such as the cc case considered in [Q. This has the consequence that the interference 
terms are numerically negative in both cases, as shown in ||: wc thank P. Zerwas for discussions on this point. 

11 There is no indication of this in the data available so far, but we do note that there are indications of other 
flavour asymmetries in the qq sea p3fl. 
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A' — 4 

A 132 — 


mi x = 200 GeV, m t ~ 2 very heavy 






squared amplitude (pb) 


negative interference (pb) 


n n 

U.U 


U.UoZ 


0.415 


-0.333 


0.2 


0.063 


0.383 


-0.320 


0.4 


0.013 


0.293 


-0.280 


0.6 


-0.043 


0.170 


-0.213 


0.8 


-0.066 


0.054 


-0.120 





A'132 = 0.4 


, m h = 200 GeV, m h = 


= 205 GeV 


\sin(f)\ 


cross section (pb) 


squared amplitude (pb) 


negative interference (pb) 


0.0 


0.082 


0.415 


-0.333 


0.2 


0.066 


0.399 


-0.333 


0.4 


0.025 


0.356 


-0.331 


0.6 


-0.016 


0.313 


-0.329 


0.8 


-0.020 


0.306 


-0.326 



Table 1: The contribution additional to the Standard Model cross section for e + e~ — > ss at Eqm = 
192 GeV, for different values of the stop mixing angle <f>. The second column gives the net additional 
contribution to the cross section due to the new terms. The third column presents the positive 
contributions due to the new squared amplitude terms in the first line of @, and the fourth column 
presents the negative terms due to the destructive interference in the second line of (|9|). 



HERA. Both the "down stop" and "strange stop" mechanisms involve the component of the 
e + beam with right-handed longitudinal polarization. This means that dramatic effects should 
be seen once longitudinal lepton beam polarization is commissioned at HERA: both models 
predict that cr(e^jo) <<< cr(e^p) and cr(e^p) <<< cr(e^p). Polarized proton beams |24| would 



also be very interesting. In the case of "down stop" model, since most models predict that 
valence down quarks should mainly be polarized opposite to the proton spin, one would expect 
that a(e^p R ) « a(e^p L )- 

The predictions for polarized proton beams in the "strange stop" model are currently less 
clear, since there is no consensus on the amount of polarization As of the ss sea in the proton. 
The most straightforward interpretation of the available data on polarized deep-inelastic £N 



scattering is that As < P5] , though this suggestion is subject to possible gluonic reinterpre- 
tation . However, the available data give no clear indication on the degree of polarization 
As/(s + s), particularly not at the large x values probed by the HERA experiments, nor do 
we know how As might be shared between the s and s in the proton. One particularly naive 
model suggests that the s quarks may have a large degree of polarization, particularly at large 



x, and that the same might also be true of the s [27[. In this case, one might expect that 



^rPr) < cr{e^p L ) , <r(e L p L ) < (?)a(e L p R ) (13) 
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with much smaller cross sections for the other beam polarizations, as mentioned above. 

It is very likely that the mystery of the large-Q 2 HERA anomaly will be resolved by the 
time such polarized experiments become possible. However, these examples serve to illustrate 
that both beam polarizations could be useful analysis tools. 
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